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ille vast majority of known non-accreting neutron stars (NSs) are rotation- 

00 

O powered radio and/or 7-ray pulsars. So far, their multiwavelength spectra 

CN 

have all been described satisfactorily by thermal and non-thermal continuum 

# > 

models, with no spectral lines. Spectral features have, however, been found in a 
handful of exotic NSs and thought to be a manifestation of their unique traits. 
Here we report the detection of absorption features in the X-ray spectrum of 
an ordinary rotation-powered radio pulsar, J1740+1000. Our findings bridge 
the gap between the spectra of pulsars and other, more exotic, NSs, suggesting 
that the features are more common in the NS spectra than they have been 
thought so far. 
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The study of neutron stars (NSs) offers a unique opportunity to understand physical pro- 
cesses under the most extreme conditions in the Universe. These super-dense, compact stars 
have extremely strong magnetic fields reaching up to a few hundred Teragauss. Thanks to the 
strength of the magnetic field and rapid rotation, young NSs are capable of generating and 
maintaining huge electric potential drops (up to a few Petavolts) across parts of their magne- 
tospheres, leading to efficient particle acceleration, pair cascades and high-energy radiation. 
The vast majority (> 98%) of known NSs are rotation-powered radio or 7-ray pulsars (1) (Fig. 
1). The continuum emission of these objects (both thermal emission from the NS surface and 
non-thermal from the magneto sphere) has been extensively studied and has provided important 
information about the structure and properties of pulsar magnetospheres, the composition and 
temperature of the surface layers, and the state of the NS interior (2). 

Recently, absorption features have been detected in the thermal X-ray spectra of a few 
unusual NSs - 'X-ray dim Isolated NSs' (XDINSs), which do not emit radio pulsations and 
likely have empty (or particle- starved) magnetospheres, Central Compact Objects (CCOs) in 
supernova remnants, which are radio-quiet NSs with low magnetic fields, and 'Rotating Ra- 
dio Transients' (RRATs), which belong to a recently discovered class of NSs seen in the radio 
only intermittently. All of these NSs are atypical compared to ordinary, rotation-powered pul- 
sars, and the spectral features have been attributed to peculiarities of these objects and exotic 
physics. For example, in XDINSs the surface magnetic field is strong enough for QED effects 
to be important, while for the very weakly magnetized CCOs the electron cyclotron energy for 
matter at the surface is low enough that it falls into the soft X-ray band. For ordinary pulsars 
(with magnetic fields ~ 10 12 G) neither of these conditions is expected to be satisfied at the NS 
surface and hence no analogous spectral features are expected in soft X-rays. 

The 100 kyr old radio pulsar J1740+1000 was detected in 2002 (3). The pulsar's spin period, 
P = 154 ms, and its derivative, P = 2.15 x 1CT 14 s s _1 , place it among other relatively young 
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pulsars in the P-P diagram (Fig. 1), suggesting that it is a fairly unremarkable pulsar, which 
loses rotational energy at a rate E = 2.3 x 10 35 erg s _1 and has a surface magnetic field 
Bq = 1.8 x 10 12 G. With the dispersion measure DM=24 pc cm -3 and Galactic latitude of 
20.3°, the derived distance of PSR J1740+1000 is 1.4 kpc or 1.2 kpc (according to the Galactic 
electron density models of (4) and (5), respectively). The off-plane location leads to fortuitously 
low interstellar absorption, making this pulsar one of the few whose spectra can be studied 
in softer X-rays (< 1 keV). For 10% ISM ionization, the measured DM implies a hydrogen 
column density N a ~ 7.4 x 10 20 cm~ 2 . The total Galactic HI column in this direction is 
« (8.2 - 8.7) x 10 20 cm- 2 (6,7). 

X-ray emission from PSR J1740+1000 was first detected in a 20 ks Chandra ACIS observa- 
tion (8). Here we report the results of observations with the XMM-Newton X-ray observatory 
European Photon Imaging Camera (EPIC) and Chandra X-ray Observatory Advanced CCD 
Imaging Spectrograph (see the supplementary text for observation details). 

The phase-integrated spectrum of J 1740+ 1000 (Fig. 2, left panel), measured with EPIC and 
ACIS, is well-fit (reduced % 2 =0.92, for v =81 d.o.f.) in 0.2-10 keV by a model containing two 
blackbody (BB) components (9) (with T£° = 71±| eV, T 2 °° = UBt\l eV and RT = km, 
= 0.647Jq 2 7 km, for d = 1.4 kpc) plus a power-law (PL) component with photon index 
T = 1.6 ± 0.6 and iV H = 9.77ji;| x 10 20 cm -2 (10) (here and below uncertainties are given at 
68% confidence for a single parameter). The isotropic bolometric luminosity of dominant cooler 
thermal emission from most of the NS surface is 2 x 10 32 erg s _1 . The hotter thermal component 
with a factor of 140 smaller emitting area could be attributed to a hot spot on the NS surface. 
Being faint and hard, the PL component contributes significantly only above 1.3 keV (Fig. 2). 
To simplify the analysis of the thermal emission, we only considered energies <1.1 keV and 
therefore ignore the PL component. The emission shows substantial pulsations with the pulsed 
fraction of 30% ± 4% in 0.2-1.1 keV (Fig. 3, left). Overall, the properties resemble those of 
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other nearby middle-aged pulsars (Geminga, B0656+14, and B 1055— 52 (11, 12)). 

Figure 3 (right) shows spectra corresponding to three selected phase intervals: [0.0-0.3] 
("dip"), [0.3-0.6,0.8-1.0] ("rise and fall") and [0.6-0.8] ("peak"). Large systematic differences 
appear between the three spectra. While overlapping at 0.3-0.5 keV and above 1 keV, the spec- 
tra differ markedly in 0.5-0.8 keV and below 0.3 keV. The relative differences in phase-resolved 
spectra cannot result from poor calibration (as all photons are detected with the same instrument 
during the same observation), and the residuals strongly suggest two phase-dependent spectral 
features. Indeed, when we added a Gaussian absorption line (GABS) to a single BB contin- 
uum (BBxGABS, where GABS=exp(-re-°- 5[(jB - jBc)/al2 ); see Table 1 caption), the differences 
between the phase-resolved spectra in the 0.5-0.8 keV range could be accounted for solely by 
variations in the GABS parameters with the BB component held fixed. The 0.5-0.8 keV line- 
like feature is strong and present in all three phase intervals, but both the line's central energy 
and strength vary with phase, becoming progressively stronger from the pulse maximum to the 
pulse minimum (Fig. 3). In the BB xGABS fits there are no systematic residuals in the 0.5-0.8 
keV range. The residuals below 0.3 keV vary in concert with the 0.5-0.8 keV absorption line 
strengths; however, with the relatively small number of photons collected we could not reliably 
characterize the parameters of the low-energy absorption feature which is at the edge of the 
EPIC sensitivity range. 

As an alternative to the absorption feature scenario, we considered a pure continuum, BB+BB, 
model. We conclude that the absorption feature interpretation is strongly preferred over the hot- 
spot scenario (see supplementary material). 

The more significant of the two spectral features has the central energy, E = 550-650 eV 
(see Table 1, bottom section), close to the crossover energy for the hot and cold BB components 
in the pure continuum model. Therefore, because even a two-component model is likely to be 
an oversimplification (one can expect the NS surface temperature to be a smooth function of the 
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magnetic co-lattitude), one could speculate that the remaining residuals could be reduced by 
increasing the level of complexity and introducing yet another BB component. However, this 
would still not account for the phase-dependent changes below ^ 0.3 keV. One could introduce 
yet another even cooler BB, but this emission would have to come from the bulk of the NS 
surface and hence should not depend on the rotation phase. Therefore, with the data in hand, 
the Occam's razor strongly favors the spectral feature interpretation. 

So far, absorption features were seen only in the spectra of atypical, exotic isolated NSs. 
These features are expected to provide an important diagnostic of the NS equation of state, 
atmospheric composition and magnetic field strength (13, 14). Various interpretations have been 
put forward. Among widely discussed options are (a) atomic transitions in the magnetized NS 
atmospheres or condensed surface layers (15, 16) and (b) cyclotron absorption features produced 
either by electrons in the NS magnetosphere (14, 17) or protons (or light ions) in a hot ionized 
layer near the NS surface (18). 

Can the explanations used for the absorption features of exotic NSs be applied to the case 
of J1740+1000? The energy of the feature we fitted is w 0.5 — 0.7 keV, while the second, 
less significant feature is at 0.1 — 0.25 keV. These energies are too large for transitions in a 
magnetized hydrogen atom in the B 1.8 x 10 12 G surface field (19). The energies possibly 
could match a transition to an excited state in hydrogen-like ions (2) for Z > 2. In this case, 
the detected spectral features would provide a direct probe of the composition and physical 
conditions of the surface layers (see e.g., (20) and references therein). 

Alternatively, the two features could be the fundamental and the first harmonic of the elec- 
tron cyclotron energy E cc = 11.6(-B/10 12 G) keV. However, in this case the absorbing/scattering 
particles must be located in the magnetosphere (21), at several stellar radii above the NS sur- 
face in a much weaker magnetic field, B(r) = B c<l (l + 3 cos 2 9 B ) 1 ^ 2 (R/r) 3 , where R is the NS 
radius and 9b is the magnetic co-latitude. Such a screen would then be analogous to the Earth's 
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van Allen belts. Such a scenario can also explain the phase dependence of the feature's width 
and central energy, because the strength of the magnetic field and the density of the absorbing 
(scattering) screen should depend on the angle between the magnetic dipole axis and the ob- 
server's line of sight. If the spectral features are indeed due to the cyclotron absorption in the 
magnetosphere, at least some regions of it must contain nonrelativistic (or mildly relativistic) 
electrons/positrons. The observed width of the 0.5-0.8 keV feature suggests that the mean ve- 
locity of these particles should not exceed 0.2-0.3 c. It is possible that a population of "warm" 
electrons exists in the closed line zone forming radiation belts (see (21) and (22)). Because the 
orientations of the pulsar spin and magnetic axis with respect to the observer's line of sight are 
not known, it is possible that the line of sight traverses some volume in the closed field line re- 
gion of the magnetosphere at all rotation phases, i.e., some absorption is always present. For the 
spin-down parameters of J1740+1000 equation 4 of (21) predicts optical depth of ~ 3 x 1CT 3 
at E = 600 eV, assuming the Goldreich- Julian charge density. However, it is widely accepted 
that the actual charge density in the pulsar magnetosphere must be larger due to pair produc- 
tion, and the multiplication coefficient of 200-300, required in our case, is comfortably below 
the typically discussed range of 10 3 -10 6 (e.g., (23) and references therein). Because it is not 
clear why heavy (Z > 2) elements would be exposed only in certain parts on the NS surface 
(thus leading to phase-dependence), we currently prefer a magneto spheric interpretation of the 
absorption features. 

There are three other middle-aged pulsars (Geminga, B0656+14, and B 1055-52) which are 
similarly nearby, bright, and have low interstellar absorption. The spectra of these pulsars show 
significant variability with pulse phase, which has been commonly attributed to the changing 
view of the non-uniformly heated NS surface (12). However, even with allowance for such vari- 
ability, noticeable phase-dependent residuals (24) in the spectral fits are yet to be explained. The 
spectral features found for J 1740+ 1000 suggest that the systematic residuals for other pulsars 
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may be better explained by a model with absorption feature(s). Indeed, there is strong evidence 
that the lines seen in XDINSs and/or CCOs also show strong phase dependence (25, 26). The 
dependence of absorption feature properties on pulsar characteristics can tell us about the phys- 
ical origin of the features. For example, in the cyclotron scenario, the energies of the lines are 
expected to correlate with the magnetic field strength and possibly with the size of the mag- 
netosphere. Atomic features would depend upon the pulsar surface density, temperature and 
composition (as is seen in, e.g., white dwarfs). 
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Parameter 




Phase interval 






Dip 


Rise and Fall 


Peak 




[0.0 - 0.3] 


[0.3-0.6],[0.8-1.0] 


[0.6-8.0] 




BB+BB (fixed T 12 = 71, 148 eV), x 


2 = 64.4 for 65 d.o.f. 


Rf ,km 


7.4 ±0.2 


7.6 ±0.2 


8.4 ±0.3 


Rf,km 


o-6i±8:82 


u - Dy -0.02 


o.8oi8:S 




BB+BB (Ti, T 2 are fitted), X 2 = 


53.7 for 59 d.o.f. 


Tf°, eV 


66 ±6 


67 ± 11 


34+1° 


Rf,km 
T 2 °°, eV 


Q 2+ 3 - 1 

y - z -2.0 

184+20 


8.4+Ji 
130+^ 


102t 4 f 
103^ 


.ffif ,km 


n S5+ - 14 

u - oo -0.10 




o n+0.4 

^■ u -o.5 




GABSxBB model, X 2 = 46.2 for 60 d.o.f. 


Tf°, eV 




94 ± 2 (tied) 




#f ,km 




4.3 ±0.1 (tied) 




£ c , eV 
cr,eV 


646+J 3 , 


635+^4 
161±fl 


548 ± 12 

qc;+22 


r 


0.92 ±0.14 


0.58 ±0.10 


0.80 ±0.57 


EW,eV 


234 ± 34 


192 ± 24 


54 ±20 



Table 1: Parameters (with la statistical uncertainties) of the spectral models fitted to the phase- 
resolved spectra in 0.2-1.1 keV range. In all cases the spectra from three phase bins have been 
fitted simultaneously for consistency. T°° (eV) and (km) are BB temperature and radius 
(as measured at the infinity), respectively. The radii are given for 1.4 kps distance. For the line, 
E c , a, r, and EW are the central energy, width, optical depth and the corresponding equivalent 
width, respectively. The column density Nu is kept constant in every case, at the value of 
9.7 x 10 20 cm~ 2 obtained for the phase-integrated spectrum (see text). EW is not an additional 
parameter, but is inferred from the other GABS parameters. 
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Figure 1: Period P versus period derivative P for pulsars of various types. Two parameters are 
most commonly used to characterize pulsars: the pulsation period P, and its rate of change P. 
The data come from the ATNF catalog (27). A number of physical indicators can be derived 
from P and P, namely the spin-down power, magnetic field, and spin-down age (over-plotted as 
labelled straight lines). Whereas the vast majority of ordinary, rotation-powered pulsars (black 
dots) are grouped in the middle of the Figure, exotic objects such as magnetars (red squares), 
RRATs (blue triangles), CCOs (red circles,) and XDINSs (orange stars) appear located away 
from the ordinary pulsar population. We have marked the NSs with reported spectral features 
(for which both P and P have been determined) with large black crosses. 
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Figure 2: Phase-integrated spectra of PSR J 1740+ 1000. Left: The count-rate spectra obtained 
with XMM-Newton EPIC and Chandra ACIS are shown together with the best-fit absorbed 
BB+BB+PL model. The errorbars shown in this and other figures are la for Poisson statistics. 
Different colors correspond to different instruments of the two X-ray observatories (as labeled). 
Here and in other spectral plots, the errorbars correspond to the data, and the lines correspond to 
the best-fit models. Right: The best-fit model photon spectrum, with the individual components 
shown (dotted lines). See text for details of the fitting and parameter values. The jumps in the 
spectrum are the edges due to the interstellar absorption. 
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Figure 3: Pulse profiles and phase-resolved spectra of PSR J 1740+ 1000. Left: Energy-resolved 
background- subtracted pulse profiles measured with XMM-Newton EPIC PN in five energy 
bands. Top right: Count-rate spectra constructed by selecting the photons from the three chosen 
phase intervals shown by red, green, and black in the left panel. We also show residuals with 
respect to the BB component (solid lines) of the best-fit BB x GABS model to emphasize the 
shape and locations of the lines we observed. Bottom right: BB x GABS fits and residuals for 
the three chosen phase intervals. No systematic residuals are seen between 0.5 and 0.8 keV. 
Here we have not fitted for the low-energy line, leaving systematic residuals below 0.3 keV. See 
supplementary material for details. 
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Figure 4: * 

Figure SI: Top: EPIC-PN count-rate spectra below 1.3 keV together with the best-fit BB+BB 
model for the three chosen phase intervals, leaving significant non-zero residuals below 0.3 
keV and in 0.5-0.8 keV range. Bottom: Confidence contours (68% and 90% for single 
interesting parameter; Ax 2 =1.0, and 2.7, respectively) for the cool BB component of the 
BB+BB fit to the phase resolved spectra. During the contour calculation, Nh was frozen at the 
best-fit value of 9.7 x 10 20 cm -2 , while the rest of the parameters were re-fitted at each point 

of the grid. 
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Supplementary Material 

Details about X-ray observations: The EPIC PN detector was operated in Small Window 
mode, providing a time resolution of 5.7 ms, while the MOS1 and MOS2 detectors were in 
Full Frame mode with a 2.6 s time resolution. After correcting for dead time, the useful PN 
exposure was 46 ks. The aim of this observation was to analyse the X-ray spectrum as a function 
of pulse phase (i.e., phase-resolved spectroscopy). We also acquired a 64.6 ks Chandra ACIS 
observation, with a much higher angular resolution of ~ 0.5" but a coarser time resolution 
of 3.2 s. Prior to fitting, all the spectra were grouped to have at least 30 counts per bin. The 
background contribution is completely negligible for the ACIS data while for the EPIC data it 
was 8%— 13% (in 0.2-10 keV), becoming significant only at photon energies E > 5 keV. In total 
we have analyzed 4,448 photons in 0.2-10 keV for EPIC and 397 photons in 0.5-8 keV for 
ACIS. Spectral modeling and fitting was done with XSPEC v. 12.7.0 (28). 

Producing pulse profile and phase-resolved spectra: Pulsars are NSs characterised by 
their very regular periodic variations in brightness. PSR J1740+1000 has a period P = 154 
ms (or, equivalently, v — 1/P — 6.49 Hz) which generally decreases with time. However, 
one can directly measure the period at the epoch of observation from the obtained data, given 
sufficiently strong pulsed signal. We used the well-known Z% (Rayleigh) test to search for the 
period in the vicinity of 154 ms (more precisely in 6.48-6.50 Hz range) in the EPIC PN data. 
The results show a clear and highly significant maximum {Z\ = 45.9) at v = 6.489649(2) Hz, 
which we used to calculate the phase = vt and fold the pulse profiles shown in Figure 3 (left 
panel) and to extract the phase-resolved spectra. 

Consideration of the alternative BB+BB model: Although the BB+BB model's fit quality 

is formally acceptable (Table 1), there are a number of serious deficiencies: 

• The x 2 values are worse than those for the BBxGABS model (Table 1). To check the 

chance occurrence probability for the x 2 reduction, we have carried out Monte-Carlo 
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simulations (see below) following the recommendations of (29). The observed x 2 im- 
provement is significant at 99.4% confidence even without taking into account the evi- 
dence for a second feature (below 0.3 keV) and a similar dependence of the strengths of 
the two spectral features on the pulsar phase (Figure 3, left panels). 

• Unlike the BB x GABS model, systematic residuals remain both in 0.5-0.8 keV and below 
0.3 keV (Fig. SI) despite the larger number of model parameters that are allowed to vary 
(12 in BB+BB vs. 1 1 in BB xGABS; cf. middle and bottom sections of Table 1); 

• the size of the emitting area of the cooler BB component fitted to the "peak" spectrum 
(R = 102l7o° km) substantially exceeds a conventional NS radius (Fig. S 1 , bottom panel) 
while fitting with a hydrogen atmosphere model would increase the apparent radius even 
further (30); 

• the "dip" spectrum cannot be fit satisfactorily by a BB+BB model. Even if we allow all 
parameters to vary, we find xl — 1-5 for 22 d.o.f, and the absorbing column becomes 
unreasonably large ( N H = 20^g 2 x 10 20 cm~ 2 ) and barely compatible with that from 
the phase-integrated fit. The intrinsically broad BB simply fails to adequately describe 
the narrow spectral peak (centered at about 0.64 keV) in the "dip" spectrum (Fig. SI, top 
panel). 

Statistical comparison of BB+BB and BB x GABS models: We calculated the ratio, r obs = 
LR BBxGABS /LR BB+BB , of the likelihood ratios for the BB x GABS and BB+BB fits to the ac- 
tual data. We then simulated 1000 "fake" spectra (using XSPEC "fakeit" command) based 
on the best- fit BB+BB model (null- hypothesis). These spectra mimic the expected noise fluc- 
tuations in a given exposure time (which was the same as for the real data). The same re- 
sponse and background files have been used. Each of the simulated spectra was then re-fitted 
with the BB+BB and BBxGABS model with the temperatures, radii and the line parameters 
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allowed to vary. For each of the fake spectra (i = 1, 1000) we also recorded the ratio, 
n = LRf BxGABS /LR B B+BB , of the likelihoods for the BBxGABS and BB+BB fits. We then 
calculated the number of instances with r\ > r obs and found only six such cases. Therefore, 
for the real data the BBxGABS model is statistically preferred over the BB+BB at a 99.4% 
confidence level even without taking into account the evidence for a second spectral feature 
with a factor of two lower energy and progressive strengthening of the features from the pulse 
maximum to the pulse minimum. 
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